Pattern formation of indirect excitons in coupled quantum wells 
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Using a nonlinear Schrodinger equation including short-range two-body attraction and three- 
body repulsion, we investigate the spatial distribution of indirect excitons in semiconductor coupled 
quantum wells. The results obtained can interpret the experimental phenomenon that annular 
exciton cloud first contracts then expands when the number of confined excitons is increased in 
impurity potential well, as observed by Lai et al. [Lai et al., Science 303, 503 (2004)]. In particular, 
the model reconciles the patterns of exciton rings reported by Butov et al. [Butov et al., Nature 
418, 751 (2002)]. At higher densities, the model predicts much richer patterns, which could be 
tested by future experiments. 

PACS numbers: 71.35.Lk, 71.35.-y, 73.20.Mf, 73.21.Fg. 



I. INTRODUCTION 

Recently great progress has been made on the research 
of exciton Bose-Einstein condensation (BEC) by making 
use of indirect excitons (spatially separated electron-hole 
pairs) in coupled quantum wells (CQW). The advantage 
of the indirect excitons is that they have a long lifetime 
and a high cooling rate. With these merits, Butov et al. 
have successfully cooled the trapped excitons to the or- 
der of IK. [l| Although there are not enough evidences 
to prove that the excitons are in BEC state, it is inter- 
esting to observe several interesting photoluminescence 
(PL) patterns. Q Independently Snoke et al. have also 
observed somewhat analogous PL patterns in a similar 
system. [3] The key features in the experiment of Ref. |2j 
are as follows, (i) Two exciton rings are formed. When 
focused laser is used to excite the sample and the prompt 
luminescence is measured in the vicinity of the laser spot, 
a ring, called the internal ring, formed. While a second 
ring of PL appears as distant as 1 mm away from the 
source, called the external ring, (ii) The intervening re- 
gion between the internal and external rings are almost 
dark except for some localized bright spots, (hi) Periodic 
bright spots appear in the external ring. Q The bright 
spots follow the external ring either when the excitation 
spot is moved over the sample, or when the ring radius is 
varied with the excited power, (iv) The PL is eventually 
washed out when the temperature is increased. 

To explain the above features, a charge separated 
transportation mechanism was proposed which gives sat- 
isfactory results to the formation of the exciton ring and 
the dark region between the internal ring and external 
ring, and even the ring in a single quantum well. [J, |5[ 
However, the physical origin for the periodic bright spots 
in the external ring is still controversial. It is commonly 
believed that the periodic bright spots are formed due to 
some kind of instability. Levitov et al. considered that 
exciton states are highly degenerate and the instability 
comes from the stimulated scattering. Q Sugakov sug- 
gested that the instability is due to the attractive inter- 



action between the high-density excitons. [7] In addition, 
Yang et al proposed that the ring bright spots result from 
the interplay between the random potential and the non- 
linear repulsive interaction of the condensed excitons. [8| 

More recently, great attention has been attracted by a 
surprising observation of trapped excitons in an impurity 
potential well. In the experiment of Lai et al. [9(, they 
used a defocused laser to excite CQW in a large area. 
Excited electrons and holes are collected by an impurity 
potential well. The area of the defocused laser spot is 
typically larger than the area of impurity potential well. 
The defoused laser spot can be applied either away from 
the impurity potential well or directly on it. It is found 
that the PL pattern is much more concentrated than a 
Gaussian with a central intensity dip, exhibiting an an- 
nular shape with a darker central region. In particular, 
with increasing the laser excitation power, exciton cloud 
first contracts and then expands. Even more interest- 
ingly, pumping by higher-energy laser, the dip can turn 
into a tip at the center of the annular cloud. In fact, this 
kind of annular shape pattern can also be found in some 
localized bright spots between inner ring and external 
ring of Butov's experiment. [2| 

Neither the mechanism of stimulated scattering [6( nor 
the mechanism of pure attractive interaction between ex- 
citons [7| is able to explain the above remarkable phe- 
nomenon. With the mechanism of the stimulated scatter- 
ing, the stimulated scattering rate should enhance when 
the particle density is increased. Consequently, the an- 
nular ring should contract all the way and no expansion 
is expected. While the mechanism of the pure attrac- 
tive interaction between excitons also has the difficulty 
in explaining the expansion of excitons. 

In this paper, we propose an alternative model to un- 
derstand the formation of the bright spots in the exci- 
ton ring. We consider that the exciton system can be 
described by a nonlinear Schrodinger equation which in- 
cludes an attractive two-body and a repulsive three-body 
interactions. The interplay between these interactions 
and the kinetic energy can lead to complex patterns, 



which are shown to be in great likeness to the experi- 
ments. The nonlinear Schrodinger equation is numeri- 
cally solved and the corresponding spatial distribution of 
excitons not only can explain the ring bright spots but 
also can describe the contraction and the expansion phe- 
nomenon of the exciton ring. Some detailed features in- 
cluding a tip occurring at the center of the annular cloud 
can also be reproduced. Our model predicts some new 
patterns at higher-density regime, which can be tested 
by future experiments. 

It is important to emphasize a key element which sup- 
ports the proposed model. As pointed out in Refs. [J, [a], 
when electrons and holes are excited by laser, they are 
hot electrons and hot holes initially. Since the drift speed 
of hot electrons are larger than that of hot holes (electron 
has a smaller effective mass), electron and holes are in- 
deed charge-separated. No true exciton is formed at this 
stage. Because hot electrons and holes have a small re- 
combination rate, they can travel a long distance from the 
laser spot. After a long-distance travel, hot electrons and 
holes collide with lattice and are cooled down. The cool- 
ing speed of electrons is faster than that of holes, and con- 
sequently cooled electrons and cool holes will meet. This 
kind of charge separated transportation mechanism has 
been used to interpret the exciton ring formation. 0, Q 
In the experiment by Butov et al. [2J, they meet in the 
region of the (external) ring. While in the experiment 
by Lai et al. 9], they meet in the impurity potential 
well. It is believed that cooled electrons and holes will 
form excitons. Thus the elemental particles in the ex- 
ternal ring [2j and in the impurity potential well [dj are 
excitons. At this stage, charges are not separated, but 
coupled or bound together. They have long life. The 
interaction between excitons results in the nonhomoge- 
neous density distribution which in turn results in the 
complex PL patterns. Since particles can only move in 
CQW, their movements are basically two-dimensional. 
PL intensity is directly proportional to the exciton num- 
ber. In the following discussion, we simply take exciton 
probability density distribution as PL distribution. 



II. REMARKS ON INTERACTIONS AND 
EQUILIBRIUM 

There are two experimental facts which are important 
to the understanding of the pattern formation of the in- 
direct excitons. (1) The distribution of the excitons is 
inhomogeneous and (2) the inhomogeneous distribution 
varies with the exciton number density. To understand 
these phenomena, the key may lie in the interactions be- 
tween excitons. First, it is very clear that the interaction 
between the indirect excitons is neither purely attrac- 
tive, nor purely repulsive. If the interaction between the 
indirect excitons is purely repulsive, it will drive the ex- 
citons towards homogeneous distribution and the exciton 
cloud will expand with the increase of the exciton num- 
ber. At present there is no experimental signature to 



show this. On the other hand, if the interaction is purely 
attractive, the system is expected to collapse when the 
exciton density is greater than a critical value to which 
there is no enough kinetic energy to stabilize the exciton 
cloud. Experimentally the collapse of an exciton cloud 
has never been observed. In addition, the case of a purely 
repulsive or a purely attractive interaction is also against 
the experimental fact that the exciton cloud contracts 
first and expands later when the laser power is increased. 
The existence of the attractive interaction does not mean 
that the exciton state is unstable against the formation 
of metallic electron-hole droplet because the repulsive in- 
teraction may dominate over the attractive one in that 
regime. 

Some remarks are in order on the existence of the at- 
tractive interaction. Different to the direct excitons in a 
bulk material or in a single quantum well, the indirect 
excitons have same polarization direction since the in- 
direct excitons are formed by electrons and holes which 
are spatially separated in different quantum wells. They 
are aligned dipoles. The interaction between the indirect 
excitons contain the dipole-dipole term and the van der 
Waals term. Q The van der Waals attraction is given ex- 
plicitly by the form -C 6 /r 6 - C 8 /r s - ■ ■ ■. When the 
spacing of the indirect excitons is large, the dipole-dipole 
interaction dominates, so the interaction is effectively re- 
pulsive. However, when the spacing of the indirect ex- 
citons becomes small, the van der Waals attraction will 
dominate the dipole-dipole interaction. It is found that 
the interaction becomes effectively attractive when the 
separation between two excitons is about 3 to 6 exciton 
radii. 7] In the current experiment, the exciton density is 
about 10 10 / cm 2 . For this density, the average distance 
between the indirect excitons is about 100 nm and the 
exciton Bohr radius as is about 10 ~ 50 nm [9( or av- 
erage distance between excitons is about 2 ~ 10 exciton 
radii. In such a case, it is reasonable to assume that the 
two-body interaction is in the attractive regime. In addi- 
tion, when two indirect excitons approach to each other, 
the exchange interaction between electrons becomes im- 
portant, which also leads to an attractive interaction. In 
fact, the attractive interaction between the excitons has 
been considered as a possible candidate to describe the 
pattern formation observed by experiments. [71 HO] 

In the dilute limit, it is reasonable to assume that the 
free-energy density of the system is given by 



F = F + V ex n - gi n 2 



g 2 n 



(1) 



in terms of the expansion of exciton density n. Here Fq 
is the free-energy density for hot electrons and hot holes 
(i.e., the case without true exciton formed), V ex is an 
external potential, and g\ and 172 are (positive) coupling 
constants associated with two-body and three-body in- 
teractions. The — (+) sign with g\ (52) term denotes 
the attractive (repulsive) nature of the two-body (three- 
body) interaction. As elaborated above, the two-body 
interaction is believed to be effectively attractive in the 
system. Nevertheless inclusion of two-body interaction 



term only is not possible to give the system a better de- 
scription (see later) . If one keeps three-body terms in (JTJ) 
as well, then it is necessary that the three-body interac- 
tion must be repulsive in order to keep the system stable. 
It will be shown soon that the inhomogeneous distribu- 
tion of the exciton can be due to the competition between 
the two-body attraction and the three-body repulsion. 

The interaction between excitons may be even more 
complex. It may include the coulomb interaction (dipolc- 
dipole interaction, van der Waals interaction) and ex- 
change interaction. It may also include electron-phonon 
and coulomb screening effect. In fact, we use two param- 
eters gi and gi to describe all the above effects. 

Another important issue is wether the excitons are in 
thermal equilibrium. If the excitons are all in the ground 
state, or a complete exciton BEC has been reached, their 
distribution must be Gaussian-like. The complex pat- 
tern observed in experiments indicates that this is not 
the case - quite a large portion of excitons are in fact 
in the excited states. Thus it is important to have 
a better knowledge on the energy distributions of the 
trapped excitons. The energy distribution in turn in- 
volves the (complex) energy relaxation and recombina- 
tion processes, which have been studied by several ex- 
perimental and theoretical groups. [11| For a relaxation 
process, when the exciton density is low (n <§C a 2 B ), the 
effects clue to the exciton-exciton and the exciton-carrier 
scattering can be neglected. In this case, the relaxation 
time is mainly determined by the scattering of excitons 
off acoustic phonons.fl2j In particular, at low bath tem- 
peratures (Tb < 1 K), this kind of relaxation rate de- 
creases dramatically due to the so-called "phonon bottle- 
neck" effects. [ll| For the recombination process, because 
the excitons in the lowest self-trapped level are in a quan- 
tum degenerate state, they are dominated by the stim- 
ulated scattering when the occupation number is more 
than a critical value. Strong enhancement of the exciton 
scattering rate has been observed in the resonantly ex- 
cited time-resolved PL experiment [I3]. Therefore, even 
though the phonon scattering rate is still larger than the 
radiative recombination rate, thermal equilibrium of the 
system may not be reached. Essentially the distribution 
may deviate from the usual Bose one.[lJ] 



III. THE MODEL 

To proceed with the above analysis, we use the fol- 
lowing effective many-body Hamiltonian to describe the 
exciton system 



H 



drip 1 * 



h 2 V 2 

2m* 



+ V ex (r) 



1> 



Si 



?± I dr^^i/jtp + — y l drip W 



W»l>, (2) 



where V^( r ) [^(r)] denotes the creation (annihilation) of 
an exciton at the position r, to* is the effective mass of 



the exciton, and V ex is the static external potential. g[ is 
the coupling constant of two-body attraction, while g 2 is 
the coupling constant of three-body repulsion. In writing 
down the above Hamiltonian, the interactions between 
the excitons are assumed to be the contacted ones (i.e., 
the s-wave approximation is assumed) . Under the mean- 
field approximation and neglecting the exciton-pair fields 
like {ipip) and (ip^ip^) and the fields involving any three 
operators, the mean-field Hamiltonian can be written as 



H 



drip' 



h 2 V< 

2m* 



V ei 



gin + g2n 



4>, (3) 



where n — n(r) = (ip'i/)) is the local density of excitons 
at r, gx = 2g' 1 , and g 2 = 9g 2 - The corresponding static 
nonlinear Schrodinger equation reads 



2to 



r VVj + {Vex - gin + g 2 n 2 )^j = Ejtpj 



(4) 



where ipj and Ej are the j-th eigenstate and eigenvalue, 
respectively. It is assumed that the system is in a quasi- 
equilibrium state, and the spatial distribution of excitons 
is given by 



ra(r) 



M 



r^OlV-iMf 



(5) 



where Af denotes the total number of energy states that 
are trapped and Wj is an appropriate probability func- 
tion associated with the energy level Ej. As a further 
assumption, we take 



<lj 



e-f 3 ^ 



(6) 



which has the form of Boltzmann distribution. Here /3 is a 
parameter used to describe the exciton distribution. The 
distribution ([6]) plays a central role in the calculation, 
which is shown to lead to qualitatively good results in 
agreement with experiments. Other possibility has been 
tested, but none of them work. 

Numerically it is convenient to first do the following 



scaling: tpj(r)/VN -v ipj(r), Ng 1 -v g 1 , and N 2 g 2 -> g-i, 
such that Eq. (|4]) remains the same look. In this case, 
n(r) becomes the probability density which satisfies the 
normalization condition J„n(r)dS = 1. Next, rescale 
tpj(r)aph — ► "0j( r ) an d r /cpL — ► r, Eq. then reduces 
to 

~ -^^^j + ( v ex - ain + a 2 n 2 )ip j = Ejipj, (7) 

where v ex = V ex /e, a x = g\j (af L e), a 2 = g 2 /(cr^ L e), 
and £j = Ej / e. Here e = h 2 /m*o~p L with ctpl being cho- 
sen by the root-mean-square radius of the exciton cloud 
observed by photoluminescence. With the above scaling, 
it is found that a 2 /a 2 = ffi/<72£ = 9i rr i*o'p lj /g2h 2 , which 
is a constant for a particular sample. In the following, we 



shall use two different values of af ja<z for the experiments 
by Lai et al. [9] and by Butov et al. [2|]. 

In connection with real experiments, three important 
points should be clarified, (i) The exciton patterns are 
fully determined by its self-trapped interaction. The ex- 
ternal potential V ex is not the main cause for complex 
exciton patterns. In the experiment by Lai et al. [9j, 
the role of impurity potential well is to collect the hot 
particles (electrons and holes) that form excitons. Thus 
we include a parabolic potential v ex for calculations in 
regards to this experiment. However, in the experiment 
by Butov et al. [2j , the ring distribution of exciton is due 
to a charge separated transportation mechanism. It is 
believed that the size of exciton ring is much larger than 
that of the impurity potential. The role of impurity po- 
tential is just leading to inhomogenence. So we will set 
v ex = for calculations in regards to this experiment, (ii) 
When an electron and a hole form an exction, it is be- 
lieved that the kinetic energy is very low. It means that 
all the excitons are self-trapped in their self-trapped po- 
tential. The particle with energy over self-trapped poten- 
tial enegy is not in the self-trapped well, and thus should 
not take this kinds of particles into account, (iii) We 
consider the excitons distributed initially as a Gaussian. 
All self-trapped eigenstates tpj (i.e. Ej < 0) along with 
n(x, y) [via Eq. (|5j)] are then calculated self-consistcntly. 



IV. RESULTS AND DISCUSSIONS 

A. Excitons in impurity potential well 

Figure [1] shows the local probability density distribu- 
tion n{x 1 y) of the self-trapped excitons in an impurity 
potential well for four different particle number a\ = 15, 
25, 40, and 55. It is natural to quote the particle number 
N by a\ since a\ oc g\ with g\ being rescaled to be ex N. 
The impurity potential is taken to be parabolic 



^ex(r) = 



-0.8 (r/CTp L - l) z 




for r < erpL, 
otherwise, 



(8) 



with a cutoff in simulating the real system. The plots are 
on a 2D plane and the color scale denotes the relative am- 
plitude of the local density n(x, y). When the irradiating 
laser power is low (therefore, the self-trapped particle 
number is few), the dilute exciton cloud is diffused be- 
cause the attractive interaction (171) is weak. When the 
laser power is increased (therefore, g± increases), stronger 
and stronger attraction drives the exciton cloud to shrink. 
When the particle number is further increased, the repul- 
sive interaction (52) becomes more important, and even- 
tually dominates over the attractive interaction. As a 
consequence, the exciton cloud expands again. The evo- 
lution of the exciton cloud with the exciton number (the 
laser power) is in good agreement with the experimental 
observation by Lai et al. [9j. 







FIG. 1: Formation of the annular distribution of the self- 
trapped excitons as a function of exciton number (a) ai = 15, 
(b) 01 = 25, (c) 01 = 40, (d) 01 = 55 with (3 = 0.0001 in 
units of 1/e and ai — 0.005af . The results are intended to be 
compared with the experiment by Lai et al. |g]. 



For a better illustration, the corresponding wave func- 
tions of the discrete energy levels related to Fig. QJd) are 
plotted in Fig. [2] The wave functions related the distribu- 
tions in Fig.QJa)-(c) are similar to those in Fig. [5] When 
the particle number is few [Fig. [1(a)], the number of the 
self-trapped energy eigenstates involved is fewer. The 
probability density distribution n is mainly determined 
by the superposition of the ground state and the first ex- 
cited states. The ground-state wave function is s-wave 
with a peak at the center [see wave function in Fig. [2(a)], 
while the first excited states are nearly twofold degener- 
ate p wave with a node at the center [see wave functions 
in Fig. [2](b) and (c)]. Superposition of these two states 
then lead to an annular distribution with a dip in the 
center. With increasing of the particle number [Fig. [T(b) 
and (c)], the number of the self-trapped energy eigen- 
states involved increases and the second excited states 
start to intervene. Second excited states are twofold de- 
generate d waves [see Fig. [3(d) and (e)]. Thus a ring with 
much higher contrast is obtained. Further increasing the 
particle number, the third excited state [see Fig. [2(f)] 
then starts to intervene the system. The superposition 
results an annual distribution with a tip at the center. 

It is useful to estimate the values of g\ and 52 with 
respect to the real system. Taking Fig. [1(b) as an exam- 
ple, it is estimated that n « 3.0 x 10 10 cm~ 2 . Since the 
experiment gives the root-mean-square radius ctpl = 10 
fj,m, the trapped exciton number N = 7rcrp L n w 9.4 x 10 4 . 
One then obtains 
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the charge separated mechanism. [3, Q The initial distri- 
bution of excitons is assumed to be homogenous in the 
ring. Considering that ctpl ~ 50/im, which is about 5 
times larger than that in Fig. [TJ so the parameter a\ is 
approximately 25 times larger than the impurity poten- 
tial well case. Figs.|3ja) -|3tc) show the formation of the 
bright ring spots for different exciton numbers and the 
ring radii (the exciton density was kept as a constant). 
With increasing the exciton number, the ring radii (so 
as the number of the spots) increases, but the density of 
bright spot remains unchanged. The periodic spots and 
the change of the spots number with the exciton number 
(or the laser power) are in qualitative consistence with 
the experimental observations. [2| In a real sample, non- 
homogeneity and some impurity may exist in the system. 
It is believed that the distorted and nonhomogeneous 
patterns result due to the impurity potential. To give 
a better fitting to the experiments, the external poten- 
tial V ex should be restored in the Schrodinger equation 
©• 




FIG. 2: The corresponding lower-energy wave functions z 
ij)i(x,y) related to the case of Fig. Hid). 



FIG. 3: The formation of the bright ring spots for different 
exciton numbers: (a) a\ — 80, (b) a± = 150, and (c) ax = 250. 
P — 10 -5 and a,2 — 10 _5 ai. The results are intended to be 
compared with the experiment by Butov et al. |2|. 



and 
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6. 54 x l(T 36 nieV. 



(10) 



Besides, if one assumes that the attractive interaction 
comes from an s-wave scattering, then the s-wave scat- 
tering length 



AttN 



2. 1 x 10 4 nm. 



(11) 



Moreover, since the exciton Bohr radius gb = 10 ~ 50 
nm, the ratio a/as — 400 ~ 2000. If we assume (3 = 
l/fceT, the exciton temperature 



According to our picture, the physical origin of the 
bright ring spots is the consequence of the competition 
between the two-body attractive and three-body repul- 
sive interactions and the kinetic energy. If excitons are 
uniformly distributed on the external ring initially, the 
attractive interaction will drive the excitons to move to- 
gether. When the local density reaches one certain value, 
kinetic energy will drive the high-density excitons to dif- 
fuse. At the same time, the repulsive interaction also 
hinders further increasing of the exciton density. The 
competitive consequence leads to form a series of clus- 
ters on the external ring. The size of these cluster is 
determined by the ratio of these three effects. 



T = 



TO*tTp L fcs/3 



0.2K. 



B. Fragmented exciton ring 



(12) 



The same physical picture can also be employed to ex- 
plain the ring bright spots observed by Butov et al. [2j 
In the present case, the exciton ring has formed with 



C. High-density patterns 

The above calculations indicate that the phenomeno- 
logical Schrodinger equation (0| gives a good description 
to the exciton distribution, both in an impurity poten- 
tial well and in the external ring. One open question is 
what will happen when the particle number in an impu- 
rity potential well is much increased. Fig. Q] shows the 
numerical results for various exciton densities. When the 




(d) (e) (f) 




FIG. 4: The density distribution for different exciton num- 
bers: 01 = (a) 20, (b) 30, (c) 60, (d) 80, (e) 200, (f) 300. 
Other parameters are the same as those in Fig. [2] 



exciton density is low [Fig. H^a)], an annular distribution 
is observed. With increasing the density of the excitons, 
the annular exciton cloud expands and a tip emerges at 
its center [Figs. HJb) & (c)]. This phenomena has been 
observed in the experiment by Lai et a/..[9J When the 
density is increased further, the annular exciton cloud 
is unstable and a series of clusters form [Fig. 0Jd)] . In 
the case that the density is high enough, the competition 
between different interactions can in fact lead to much 
complex patterns beyond the present experimental ob- 
servations. The clusters eventually form a pattern of tri- 
angle lattice [Figs, life) & (f)]. We point out that all these 
patterns could be observable under current experimental 
conditions when the density is large enough 



low while luminous efficiency is high, excitons can not 
reach the thermal equilibrium state. Due to the com- 
petition between the self-trapped and kinetic energies, 
complex exciton patterns occur (as discussed above). 
When the temperature is increased, excitons are not fully 
cooled and correspondingly self-trapped interaction con- 
fines only part of the excitons. The attractive interaction 
can not compensate the exciton kinetic energy and exci- 
tons will distribute homogenously in 2D plane. In this 
case, the pattern is washed out. If the temperature is 
higher than the indirect exciton binding energy ~ 3.5 
meV, [la. Il6| most of excitons become ionized and are in 
a plasma state. No pattern can be observed in this case. 
In order to realize exciton BEC, the further experimen- 
tal work should look for effective method to obtain the 
excitons with low combination rate and short relaxation 
time. The further theoretical work should focused on the 
exciton interaction beyond the mean field approximation. 

To summarize, we have demonstrated that the exciton 
distributions observed in experiments can be explained 
by the competition between the self-trapped interaction 
and the kinetic energy. A nonlinear Schrodinger equa- 
tion including short-range two-body attractive and three- 
body repulsive interactions is used to describe the exciton 
behavior. The interplay among the two-body interac- 
tion, the three-body interaction, and the kinetic energy 
not only explains the experimental observations, but also 
leads to rich patterns, which could be tested in future ex- 
periments. 
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